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Abstract 

For multi-level time-dependent quantum systems one can construct superadiabatic 
representations in which the coupling between separated levels is exponentially small 
in the adiabatic limit. Based on results from |BeTei| for special Hamiltonians we 
explicitly determine the asymptotic behavior of the exponentially small coupling term 
for generic two-state systems with real-symmetric Hamiltonian. The superadiabatic 
coupling term takes a universal form and depends only on the location and the strength 
of the complex singularities of the adiabatic coupling function. 

As shown in BeTei, , first order perturbation theory in the superadiabatic represen- 
tation then allows to describe the time-development of exponentially small adiabatic 
transitions and thus to rigorously confirm Michael Berry's {Berj, predictions on the 
universal form of adiabatic transition histories. 
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1 Introduction and the main result 

We consider the dynamics of a two-state time-dependent quantum system with state vector 
■0 G described by the Schrodinger equation 

{iedt-H{t))i;{t)=0 (1) 

in the adiabatic limit e ^ 0. The Hamiltonian H{t), t € M, takes values in the real- 
symmetric traceless 2 x 2-matrices of the form 

TT(i\-l( cos6l(t) sin6l(t) \ 
^^^>- 2^ sm9{t) -cose{t) J ■ 



This is the prototype of all adiabatic problems in quantum mechanics and the "Adia- 
batic Theorem" states that if H{-) £ C^(M) the system can be decomposed into two 
scalar equations which are decoupled up to errors of order e. Under further regularity 
assumptions on H(t) the error bound can be improved and for suitable real analytic H(t) 
it is even exponentially small in i. This adiabatic decoupling is at the basis of under- 
standing a large number of physical phenomena related to the separation of time scales, 
ranging from the classical Stern-Gerlach experiment for the measurement of spin to the 
dynamics of molecules; see (PSTt .Te for recent reviews. Despite its asymptotic smallness, 
the exponentially small coupling that generically remains has itself important physical 
consequences such as the non-radiative decay in molecules. In the scattering limit the 
exponentially small non-adiabatic transitions are quantified by the Landau-Zener formula 
and its generalizations, with rigorous justification given in [jKPllJoj . In this work we treat 
the problem of explicitly determining the exponentially small non-adiabatic coupling for 
arbitrary finite times t in order to obtain a complete understanding of the nature and 
the time-development of non-adiabatic transitions. Our main result is the construction 
of a family of unitary maps [/""(t) that brings into almost diagonal form © with 
off-diagonal elements c^^it) that are exponentially small in e and are explicitly given at 
leading order. The construction works under assumptions satisfied for "generic" Hamil- 
tonians in a sense to be made precise. Our work was motivated by results of Berry |Berj 
(see also BerLi , LiBej ). He argues that the time-development of non-adiabatic transitions 
is determined solely by the complex singularities of 6'{t) closest to the real axis, which 
are "generically" first order poles. For such generic poles the transition histories then 
have the universal form of an error function. Despite substantial progress in adiabatic 
theory during the last decade, e.g. |JoPtl iNel iMa] . a rigorous justification of Berry's 
conjecture for the generic case remained an open problem until now. We are aware only of 
two results |HaJo[ |BeTei| , which both deal with special and "non-generic" Hamiltonians. 
However, the present work is a continuation of |BeTei| and relies on techniques and results 
developed there. We also refer to |BeTei| for a more detailed introduction and a guide to 
the literature on adiabatic theory in quantum mechanics. 

Before we describe our result in detail let us shortly comment on the special form of 
the Hamiltonian (jSJ, whose eigenvalues are equal to it^ independent of t. Berry and Lim 
[BerLij observed that the Schrodinger equation 

{ieds - H{s))i^{s) = (3) 

for any traceless real-symmetric Hamiltonian H{s), 

cos 9{s) sin^(s) \ , , 

sin^(s) -cos^(s) y ' ^ ^ 

with eigenvalues p±{s) = ±p(s) = ±^^/X'^(t) + Z'^{t) that satisfy p+{s) — p-{s) = 2p{s) > 
g > for all s £ M, can be brought into the form & (|2j through the invertible 
transformation to the natural time variable 

t{s) = 2 [ p{u) du . (5) 
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While our main result is formulated for and @, the difficulty of the problem stems 
in parts from the fact that our assumptions must be general enough to be satisfied by 
Hamiltonians © that arise from generic analytic Hamiltonians of the form Q through 
the transformation Q. E.g., it turns out that these assumptions prevent the use of 
standard Cauchy estimates. Note that, as was observed by Berry jBer j . Equation with 
a complex-hermitian Hamiltonian is unitarily equivalent to a similar equation with a real- 
symmetric but e-dependent Hamiltonian. Our results on & ((2)) are sufficiently uniform 
to also apply to this case and thus cover generic analytic self-adjoint 2 x 2-Hamiltonians, 

cf. IBeTeal - 

It is well known (see |BeTei| for details and references) that for all n G No there is 
basis transformation such that in the n**^ superadiabatic basis the off-diagonal elements of 
the Hamiltonian are of order e""*"^, i.e. there is a unitary U^{t) such that 



U^{t) (iedt - H{t)) U^itf U^{t)m =■■ Mt - H:it))r{t) = 



H^(t)= \l^^^. "7. h c^(t)=0(e"+^), and ^"(t) = ^^(tWt) • (6) 




with 



m 

Here p'^{t) = ^ + 0(e^)- Note that U^{t) =: Uo{t) is the orthogonal transformation that 
diagonalizes the symmetric matrix H{t). It is independent of e and maps to the adiabatic 
basis. However, in general, lim„^oo Ic^ (i)| = oo for all e > and the coupling can not 
be eliminated completely for fixed e by going to higher and higher superadiabatic bases. 
Instead, for each e > there exists an optimal = n(e) for which n i— > |c"(t)| attains it 
minimum |c"^(t)|. This defines the optimal superadiabatic basis. In order to determine 
and (t) it is necessary to understand precisely the asymptotic behavior of the coupling 
(i) as n ^ cxD. 

As in the case of Berry's non-adiabatic transition histories one expects that for large 
n the superadiabatic coupling function c"(t) is determined by the singularities of the 
adiabatic coupling function c^{t) = ^0'{t) and thus has a universal form. From the 
abstract asymptotic analysis point of view this universality is just another manifestation 
of Darboux' Principle, which in its original form says that the late coefficients in the Taylor 
series of an analytic function are determined by the convergence limiting singularities, see 
Theorem [SJ Dingle |Dij realized that the same idea also applies to various divergent series 
arising for example from asymptotic expansions of integrals. In our case the coefficients 
(t) are determined by a non-linear system of recurrence relations starting with ^0'{t) and 
involving differentiation, integration and multiplication of terms. We show that Darboux' 
Principle can be applied also to this system of recurrence relations, i.e. that the large n 
asymptotics of the coefficients c"(i) depend solely on the convergence limiting singularities 
of 9'. As to be discussed in Section 2, our system of recurrence relations can be interpreted 
as the formal asymptotic expansion of the solutions to a system of ODEs. It might well be 
that established techniques of asymptotic analysis can be used to determine the asymptotic 
behavior of c"(t) as n ^ oo by studying this system of ODEs, but this is far from obvious. 
Instead, our proof is based on a direct analysis of the recurrence relation using a family of 
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norms tailored to Darboux' Principle and introduced in Definition ^ below. One merit of 
our approach is that, in principle, it can be applied also to more complicated recurrences 
as arising, e.g., in the context of constructing precise coupling terms between different 
electronic levels in molecular dynamics, see |BeTe3| . 

We now describe our main result in detail. Since our construction is local in time, we 
can restrict our attention to a compact interval I C M. A sufficient condition for our main 
theorem is that the singularities of the analytic continuation of 9' are of the form 

N 

e'{z - zo) = + Y.{z - zor^^h,{z - zo) (7) 

where |Imzo| > 0, 7 G M, < 1 and hj is analytic in a neighborhood of for j = 1, . . . ,N. 
Following arguments of Berry and Lim |BerLij we show in Section 4 that this condition is 
fulfilled for generic Hamiltonians of the form (jlj). However, in all but a few non-generic 
cases, —aj ^ N for at least one j. The technical problem arising from this fact is that 
by removing the leading singularity one does not obtain a function which is analytic in 
a larger region. As a consequence, it is not sufficient to use standard Cauchy estimates 
to show that the remainder terms are asymptotically smaller than the contribution from 
the leading singularity. Instead we introduce the following norms tailored to Darboux' 
Principle. 

Definition 1. Let tc > 0, a > and / C M be an interval. For / G C°°(I) we define 



{I,aA) '■= sup sup 



and 

F^,tM) = {f^C^{I)--\\f\kl,a,t.)<^}- 

The connection of Definition ^ with ((7|) is given by Darboux' Principle, Theorem [51 
which allows to translate the information about the complex singularities of 6' into infor- 
mation about the late coefficients of the Taylor expansion of 6' on the real line. Taking 
/ = {t} we obtain 

T{a + k) 
t 



= C < 00 \f^'\t)\<C^^^ VfcGN. (9) 



Consequently ate) < ^ some a > implies that / is analytic at t and that 

the Taylor series at t converges at least inside the disk Dt^{t) of radius tc- Suppose that 
the Taylor series has finitely many singularities on dDt^ (t) , all of them being of the form 
(z — zo)~^^hk{z — Zo)-, hk analytic near the origin, ak > 0, then Theorem El implies 

\\f\\({t},f3A) < °° ^ ^ - maxofc . 

Remark 1. One might be tempted to think that for functions / that are analytic in 
Dt^{t) the norm ^tc) equivalent to 

Cauchy ._ „ | j! u i _M ~ \z\) 
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However, standard Cauchy estimates only yield 



which is larger than @ by a factor of + q. There may be ways to improve, but not up 
to equivalence of the norms: for the elliptic theta function ^3(2:) = Yl^=o^"' 1 obviously 
11^3 II (Oo 1) < 00 if and only if a > 1. On the other hand, an elementary estimate shows 
that ll^sll^QT'^y < 00. The reason of the discrepancy is that the Taylor coefficients of 

functions with a dense set of singularities on the boundary of the disk of convergence (as 
^3 has) have worse asymptotics than those of functions with isolated singularities. In many 
problems of asymptotic analysis this lack of preciseness of || • ||?|j|'q i ) V^^Y^ 110 role, since 
the leading singularity is isolated. Then one can subtract the leading singularity and the 
remainder is analytic on a slightly larger domain. In that case Cauchy estimates applied 
to the larger domain yield sufficiently small error terms. However, in our case the form 
of the function near the singularity requires the use of the precise norms ||/||(/ q, t^), 
since subtracting the leading singularity does not increase the domain of analyticity. 

The norms |H|(/Q,t^) have very convenient mapping properties under differentiation, 
integration and multiplication of functions, which are summarized in Proposition ^ and 
Proposition |^ These mapping properties are the key ingredient for our analysis of the 
recurrence relations defining c"(t), of which we will now give precise assumptions and 
results. For 7, 1^-, ic G let 

e'^{t) = \-f( — — 

be the sum of two complex conjugate first order poles located at tr i itc with residues =Fi7- 
Then, as to be discussed in Section 4, for zq = + itc condition generalizes to 

Assumption 1: On a compact interval I C [tj. — tc, tj. + td with tj. £ I let 

e'{t) = e'^{t) + e[{t) with e[{t) ^ F^^t^i) (ii) 

for some 7, tc, G I^; < a < 1. 

It turns out that under Assumption 1 the optimal superadiabatic basis is given as the 
n*^ superadiabatic basis where < iTg < 2 is such that 

ne = — — 1 + (Te is an even integer. (12) 

Our main result is the leading order asymptotics of c^^{t) for t £ I. For times t that do 
not belong to an interval satisfying Assumption 1 we establish bounds on {t) which are 
exponentially smaller than the exponentially small leading order terms near the singular- 
ities. For this we assume 
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Assumption 2: For a compact interval I and some t >tc let 0'{t) G Fibril)- 

The case of degenerate / = {t}, t G M, is exphcitly allowed. Assumptions 1 and 
2 are formulated in such a way that, in principle, 6' need only be known on the real 
axis. However, in practice we will check these assumptions by analyzing the complex 
singularities of the analytic continuation of 9' , cf. Section 4. 

In our main theorem we do not only control the asymptotic behavior of the Hamilto- 
nian in the optimal superadiabatic basis, but we also obtain constants which are uniform 
on compact intervals of the other parameters tc, a and 7. This makes the formulation 
somewhat involved, but is necessary, e.g., for the study of hermitian but not symmetric 
Hamiltonians and for the study of the scattering regime. 

Theorem 1. Let Jt^ C (0,oo), C (0,1) and Jy C (0,cxd) be compact intervals. 

(i) There exists Eq > and a locally bounded function (j)2 ■ with (j)2{x) = 0{x) 

as X ^ 0, such that for all H(t) as in ^ satisfying Assumption 2 with t^ G Jt^, for 
all £ G (0, eo] Oj'n-'i all t ^ I the elements of the optimal superadiabatic Hamiltonian 
^ and the unitary U^^{t) with as in M'^) satisfy 



and 



\U^^{t)-Uo{t)\\<e<P2{}\e'\,^,^,)). (14) 



(ii) Define 

cAt) = 2i 75 sm (Z^) e-^ e-^ cos - ^ + ^) • 

There exists eo > and a locally bounded function (pi : M"*" — > M"*" with (j)i{x) = 0{x) 
as X ^ 0, such that for all H(t) as in ^ satisfying Assumption 1 with tc G Jt^, 
a G Ja, 7 G J-y, for all e G (0, eo] and all t £ I 

|c^^(t)-c,(t)| <ei-"e-T0i(M), (15) 

where M = max | , ||^rll 

Remark 2. Assumption 1 implies Assumption 2 on the same interval / with t = tc- 
Hence, Assumption 1 implies also and H14|) with t = tc- Furthermore this shows 
that in this case the bound on in (|13)) is optimal with respect to the dependence on e. 
Again, this is only possible since we use the precise norms ||/||(/q,j^) instead of standard 
Cauchy estimates. 

Remark 3. Note that the explicit term in 1)15(1 is asymptotically dominant only if |t — tr| = 
0{^/£). Since typically r > tc in Assumption 2, for all other times t the bound given in 
(|13|) is asymptotically smaller than the error term in H15|). 
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Remark 4. It was shown in |BeTei| how to derive from Theorem ^ using first order 
perturbation theory in the optimal superadiabatic basis, the universal transition histories 
predicted by Berry |Berj . 



For generic analytic Hamiltonians the whole real line can be covered by intervals satis- 
fying either Assumption 1 or Assumption 2. Under additional conditions on the location 
of the singularities of 6' , we can also consider the scattering problem and recover the well 
known Landau-Zener formulas for the adiabatic transition amplitudes. Then the decay of 
the exponentially small coupling for large times can come either from or from 

the r-dependence of the exponent in H13|) . 

Our paper consist of two parts. The main part and the key mathematical point of our 
work is the proof of Theorem ^ Our proof relies on our previous results in |BeTei| , where 
we established Theorem ^ assuming 9'{t) = ^o(i). In Section 2 we recall the necessary 
tools and results from |BeTei| and prove Theorem Q postponing the proofs of the key in- 
equalities to Section 3. The main mathematical challenge is to determine the asymptotic 
behavior of the solutions of a system of recurrence relation, which, as shown in |BeTei| , 
yield the couplings c"(t). This is done in Section 3 and the analysis heavily relies on map- 
ping properties of the norms || • ||(7,a,tj,), which we believe are of independent mathematical 
interest. In Section 3 we also use a combinatorial lemma, whose rather involved proof is 
postponed to the Appendix. In Section 4 we finally discuss several issues concerned with 
the transformation (jSJ. In particular we present the argument of Berry and Lim [BerLij 
showing that Assumption 1 is "generically" satisfied. A more detailed analysis of this 
point as well as an analysis of interesting non-generic cases and of the scattering problem 
are postponed to |BeTe2| . This is because the mathematical problems involved are of a 
completely different type from the main problem solved in this paper. 

Acknowledgements: We are grateful to Vassili Gelfreich for several helpful remarks. 
We also profited from discussions with Gero Friesecke, Alain Joye and Florian Theil. V.B. 
thanks the Mathematics Institute of the University of Warwick for hospitality and the 
Symposium "Mathematics of Quantum Systems" organized by G. Friesecke for financial 
support. 

2 Superadiabatic representations and optimal truncation 

In this section we prove Theorem^ The mathematical object to control is the Hamiltonian 
© in the superadiabatic representation. This can be achieved by studying superadiabatic 
projections. In the simple model at hand, our understanding of these projections and 
their relation to the unitary is rather complete and has been described in |BeTei| . For the 
convenience of the reader, we give a synopsis here. 
The n*^ superadiabatic projection 



n 




(16) 



fc=0 
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is the unique operator (which is a 2 x 2 matrix in our case) with 

{T,(n)f _ ^(n) ^ 0(£n+l) ^ud (17) 

ieSf -i/,7r(")j = 0(e"+i) (18) 

for all 71 G N. Here, \A, B\ denotes the commutator of the matrices A and B; ttq is the 
adiabatic projection, i.e. the projection onto the eigenspaces of H. tt^ can be constructed 
recursively by using the basis 

^ = f ? "M ' Y = -2H, Z = -Y'/9', W = 1 



.1 . 
q£ ]^2x2 g^j^^ making the Ansatz 

TTfc = XkX + y^Y + ZkZ + WkW . 

It turns out that = for all k, while the remaining coefficients fulfill the recursive 
differential equations 

XI = -^e', yi = zi = (19) 

and 

Xn = -i(4-l - 0'yn-l), (20) 
ra-1 

Vn — ^ ^ ( XjXn—j ~\~ VjUn—j ~l~ ZjZn—j)^ (^1) 

i=i 

Zn = -ix'n-l- (22) 

In addition, the differential equation 

y'n = -G'Zn (23) 

holds for each n G N. 

In [BeTeil we construct a unitary matrix U^{t) which diagonalizes the self-adjoint 
matrix 7r^'^^{t) and achieves 



U^{t) {iedt-H{t))Ur{t)=iedt 

with 



p^(t) = 1/2 + 0(^2) and c^(t) =e"+i(x„+i(t)-z„+i(t))(l + 0(e)). (24) 
In Theorem 131 we will prove that under Assumption 2, 

\xn+i{t)\ , \zn+m < ll^'ll^^) (^exp(42 \\e'\Q - (25) 



8 



where H^*'!!)^!) = ||^'||(7 i ,-)• Note that the right hand side above is as 

uniformly in r > tc ^ inf Jt^ > 0. Using H25() and the corresponding inequahty for i/n+i 
from Theorem El in the exphcit formulas given in Section 3 of iBe^lteil , it is not difficult to 
see that ()14() holds, and that 1)24(1 can be sharpened: there exists a locally bounded function 
(j) with (j){x) = 0{x) as X — > 0, such that for all H{t) as in Q satisfying Assumption 2 
with tc G Jtc and for all t € / 



\p^it) - 1/2| < e^^dl 



1(1)^ 



and 



Combining (|25l) and (|27() . we obtain 



(26) 



(27) 



|c^(t)|<e"^ 



n! 

rn+l' 



l,r))' 



l(/,l,r) 



where i;^ has the same properties as (p, uniformly in the class of Hamiltonians just discussed. 
To arrive at ([T3|l. we take = tc/e, use Stirling's formula and analyze the asymptotics 
of the terms involved. The procedure is performed in detail in |BeTei| , and from the 
calculations there it is again obvious that uniformity in the Hamiltonians is not lost. The 
only trivial difference is that since we do not truncate at the optimal value r/e of n but 
rather at tc/s, we obtain in only a factor of exp ( — ^(1 + In j^)) instead of exp ( — ^)- 
Thus we have shown part (i) of Theorem ^ 

As for part (ii), let M be defined as in Theorem ^ In Theorem 0] we will show that 
there exists a locally bounded function (pi with (pi{x) = 0{x) as x — > 0, such that for all 
H{t) as in Q satisfying Assumption 1 with tc £ Jt^, Oi £ Ja, ^ J^, and for all t G / 



2sin(77r/2), 



n+l 



IT 



Re 1 + i 



t - tr 
tc 



< 



t 



n+l 



MM) 



provided 6' fulfills Assumption 1 and n is even; Zn+i = in that case. Now ((2 
optimal truncation show (|15|). and the proof of Theorem^ is finished. 



(28) 



281) and 



Remark 5. In jBeTeij , we used p3|) and converted the nonlinear recursion into the linear 
but nonlocal recursive integro-differential equation 



— Zn+2 



II I 



0' Zn ds. 



(29) 



Since we treated the special case where = in Assumption 1, the calculations were 
rather explicit and we obtained the analogue of Theorem ^ with even better error bounds. 
In the general situation, there is no way to avoid the nonlinear recursion (but even so, 
(|29|) will be useful). As an added bonus of not resorting to (|29|) . all our results are local. 
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Remark 6. As pointed out to us by Vassili Gelfreich, (|19|) - H22() is connected to the set of 
singularly perturbed algebraic-differential equations 



dtX{e,t) 
dtZ{e,t) 



ieZ{e,t), 

ieX{e,t) -6'{t)Y{e,t), 
e{-X\e,t) + Y\e,t) + Z\e,t)) 



with the initial condition X{0, t) = 9'{t),Y{0, t) = Z{0, t) = 0. Indeed, consider the formal 
series expansion for X, Y and Z, i.e. 



with similar expressions for Y and Z. Then H19|) - H22() are just the equations for the 
coefficients of the expansion. This opens the possibility to treat the problem using e.g. 
Borel summation, but it is not clear to us whether this would be successful. On the other 
hand, given the connection above, it may well be that our approach, to be presented in 
the following section, can be used successfully in the theory of singularly perturbed ODE. 

3 Solving the functional recursion 

In this section we examine the recursion (|19p -(|22 |) and prove, in particular, the estimates 
H25j) and (|28|) . The main ingredient to our proofs is the family of norms from Definition ^ 
Recall that for > 0, a > and a compact interval I we defined 



oo 



X{e,t) = ^e^xk+i 



k=0 



\\fkl,a,t.) :=supsup d'^fit) 



r{a + k) 



c 



< OO 



(30) 



and 



F^,tAl) = {f^C^iI)--\\f\\iI,a,t.)<Oo] 



Often tc and I will be fixed, and then we will simply write ||.| 
mapping properties of are crucial. 

Proposition 1. Let tc > be fixed, a,(3 > and t G M. Then 



and Fa- The following 




b) 11/ 



{/,Q + l,ic) 



<ll/ll(/,«,tc)- 



c) Let B{a,l3) = T(a)T{P)/T{a + /?) denote the Beta function. Then 
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Proof, a) and b) follow directly from the definitions. Turning to c), for A; > we have 



j2(':)r{a + l)T{P + k-l). 
1=0 ^ ^ 



We thus have to investigate the sum in the last line above and relate it to T{a + (3 + k). 
To do so, we use a nice trick, which is presumably well known. For — l/2<t<l/2 1et 



Then d^'hp = h/^^n- Now consider d''{hahp). Then on the one hand. 



l-t 



= B{a,P)ha+i3+k- 



On the other hand, of course 



Now we take t = Q and use /i/3(0) = r(/3). Then the above calculations give 
B{a, /3)r(a + [5 + k) = Y,()^^{oi + l)T{P + k-l). 



Inserting this in the calculations from the beginning of the proof of d), we find 



< 



II II r(a + /? + fc) ^^ 



for each k>0, and consequently 

\\f9\\(a+p) <^K/5)lM 11(a) liyil(/3) 



By taking ^ = 1 in c) we arrive at 



<t': 



Via) 



We will also need the following somewhat more special property of the norms: 



(31) 



□ 



(32) 
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Proposition 2. Let s e I and a > I. If f e Fa^t^I), then t f{r)dr G Fa„i,t^(I), 

and 

'(a-l)\t\ 



/(r)dr 



< max 



(a-l) 



tr 



1 



In case a > 2 and \t — s\ < tc this simplifies to 

ft 



/(r)dr 



<(a-l) 



(a-l) 



(a) ■ 



Proof. We have 



/(r) dr 



< t ■ 



< \t 



and for k > 1 



sup 

k>l 



fis)ds 



t. 



{a-l)+k 



r(a -1 + k) k>o 
The claim now fohows from the definition of 



d'f 



r(«) 

(a) ■ 



j.a+k 



oo r(a + k) 



(a) ■ 



(a-l) 



and the fact r(a) = (a — l)r(a — 1). 



□ 



Remark 7. The intuition behind the norm 



(/ a tc) ^® that when it is finite, the function 
/ behaves equally good or better than the function t i— > ^.^ ^^^^ when taking derivatives. 
The amazing and useful fact stated in Proposition ^ c) is that multiplication not only 
leaves this property intact, but even furnishes a factor that becomes small when either a 
or P become large. It is this property that gets all our estimates going. 



Theorem 2. Suppose that Assumption 2 holds, and write 

■= 11^11 (/,l,r) <°°- 

Then for each n G N, 



"Il{7,?i,r) 



< \\e\^(eM42\\o'\\l))-l) 

< 11^11(1) (exp(42||^'||^,p-i) 



|ynll(W) < ;73lH(42 11^11(1)) -1 



(33) 
(34) 
(35) 



Remark 8. The Douglas- Adams-constant M = 42 comes out of our proof in a natural 
way. Numerical calculations suggest that Theorem |21 holds with M = 1, but this is 
probably much harder to prove. There is also numerical evidence that the asymptotic 
behavior for large is not optimal. It appears that exp(M ||^'||^/^^) is still an upper 

bound, while exp(M ||^'||(i)) is not. 
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Proof of Theorem\^ We define C„ and D„ recursively through Ci = /2, Di = 

and 

!C„_i for n even, 

L-n-i + lor n odd, 

n - / Yl=iB{k,n - k){CkCn~k + DkDn-k) for n even, , . 

" \ for nodd. ^''^^ 

We now show that for each n G N, 

\\Xn\\{I,n,T) < C'n, (38) 

ll^'i|l(/,n,r) - (-^9) 

||yn||(,,„,.) < Dn. (40) 

This is checked directly for n = 1 and n = 2. Suppose it holds for n — 1. If n is even, then 
a;„ = so H38() trivially holds, and ()22() implies 

ll^"ll(n) = Ikri-lll(n) ~ II ^"-1 II (n-1) — C^n-l = Cn- 

Proposition ^c), 1)21(1 and the fact that either j;„ = or = at any given n G N yield 

n-l 

||yn||(„) < ^{\\XkXn-k\\(^n) + \\ykyn-k\\(^n) + \\^kZn-k\\{n)^ 
k=l 
n-l 

< ^B{k,n- k){CkCn-k + DkDn-k)- 

k=l 

If n is odd, it follows from 1^0]) that Vn = Zn = 0, and 



1 



< ||2;n-l||(„_i) + -^^—^ P'\\{1) l|yn-l||(„_i) < Cn-1 + '-^^^^Dn-l- 

This proves ()38j) - ()4fl() . From (|36() it follows immediately that 

l^n — L'n-l H — H ^--L^n-l H — 

n — 1 n—1 n — 2 

Thus it is sufficient to control the Dj. We claim: 
Lemma 1. For M > 42 and all even n G N, 



n—1 ^-^ V. 
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The proof of this purely combinatorial fact is somewhat involved and deferred to the 
appendix. From 1)42(1 and (|4U() it now follows immediately that 



\yn\\{n) 



< Dr, < 



exp(M||^'||2^))-l 



n — 1 



for all n G N. Using (jSHI) and (HTJ) we obtain 



\Xr, 



l(n) 



< C„ < ||^'||(,) U + (exp(M \\e'\\l^) -1)J2 ^ 



< ll^'ll^^^ (exp(M 11^' 



1(1)' 



and the same estimate applies to ||-2n||(„)- The proof is finished. 

Proposition ^ a) now immediately implies 
Theorem 3. Suppose that Assumption 2 holds and write 



< 



□ 



ir 11(1) ir ii(/,i,r) 

Then for each n G N and each t GM, we have 



SUp|Xn(t)| < 

t&i r" 

(n - IV 

sup|zn(t)| < ^ ^ 



1(1) 



1(1) 



exp(42 \\d' 



exp(42||6l' 



|2 ^ 

1(1)' 

|2 ^ 
1(1)' 



sup|y„(t)| < l!i_^(exp(42||0'||J^p-l). 



It is interesting that although the inequalities in Theorem were derived using some 
seemingly rather crude estimates, they are optimal up to constants in the two most im- 
portant asymptotic regimes: For large n and each ||^'||(-^) as well as for small and 



1(1) 



each n the results in IBeTeij are an example that displays exactly the asymptotic behavior 
predicted by Theorem |31 

Surprisingly, the accurate asymptotics of the recursion (|T^ - (P^ under Assumption 1 
are not difficult to obtain from the results of |BeTei| once we have the uniform bounds of 
Theorem 131 and use our norms. 

Theorem 4. Let Jt^ C (0, cxd), C (0,1) and C (0, oo) he compact intervals. There 
exists a locally bounded function <j)i : with (piix) = 0{x) as x — > 0, such that 

for all H(t) as in ^ satisfying Assumption 1 with tc £ Jt^, a G Ja, 7 G J-y and all t G I 
we have 



Xn{t) - ic. 



(n- 1)! 



Re 



Zn{t) + C. 



(n- 1)! 



fn 

''c 



Im 



1 -i- 



1 -i 



.t - tr 

t - tr 
tr 



< 



(n-l)\ 



-n 



-l+a 



MM), 



< l!^Jl:n"i+"</,i(M), 



fn 

''c 



(43) 
(44) 
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where M = max | , and = lERtiM, 

Proof. Let Xnfl, yn,o and Znfi be defined via the recursion (^01) - (|^^ started with xi^o = 
i0Q/2. This is the situation where 6' just consists of a pair of simple poles, and Theorem 3 
of [BeTeil implies and for Xnfl, Unfl and Znfl and any a < 1. Uniformity in 
the parameters 7 and tc is not spelled out there, but again it is easy to derive from the 
estimates given. Let us now write 

= Xn — Xnfl, Vn = Un ~ Unfl and Cn = ~ Znfl- 

The proof will be done as soon as we show 

knMCnl <^— ^n-i+>i(M) (45) 

''c 

for some (j) with the properties given in the Theorem, uniformly in the parameters a and 
tc- Without loss we assume = 0, and we write instead of Fk^td^) etc. The main step 
is 

Lemma 2. ^„ G Fn-i+a and Qn G for each n G N, and there exists (j) : 

with (j){x) = 0{x) as X — > such that 

||Cn||(„_i+„)<0(M) and ||en||(„_i+„)<</'(M) 

for all n S N, uniformly in a £ Ja, tc £ Jt^. 

Proof of Lemma\^ We first prove the assertion for the Cn- For odd n, ^„ = and there is 
nothing to prove, so let n be even. Since C2 = ^rV^' assertion is true for n = 2. Using 
(|23)) along with the recursion, we find 



Zn+2 = < + {O'fzn + 0" (^J^ e'Zn ds + y„(0) 



(46) 



which is obviously just another way to write (|29() . We decompose ()46p into terms which 
contribute to Zn+2fl and those that do not, with the result 



Zn+2 = z'nfl + (^o)^^n,0 + Oq (^j O^Zufl ds + y„,o(0)^ + 



(47) 



= ^n+2,0 

+Cn + {O'fQn + {2e'^e[ + {e[f)Znfl + (48) 



+e" f e'Cn ds + 9'^ f e[znfl ds + e'; f e'z^A ds + 0"y„(o) - ^[;y„,o(o). 

JO Jo Jo 

The terms in (|48() and (|49() contribute to Cn+2; and we are going to estimate the ||. 
norm of each of them, using Propositions ^ and |21 Starting with (|48j) . we have 



(49) 

\[n+l+a)' 



IK"ll(n+l+a) - IKII(n-l+Q)' (^0) 
|(^')'Cn||(„+i+,) < i3(2,n-l + a) 11(^01(2) IICII(„-l+.)< 

^ (n + a)(n-l + a) H^, \\^kn-...V 
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and 



\\{2e'^e[ + {e[?)zn4(n+i+a) < b{i + a,n) \\9',ie', + lkn,oll(„) 
r(a)r(n) ,„ 



< 



r(n + a + l) II "^"H Vll 011(1) 
Turning to (|i9|) . let us first note that 



+ 



(1)7 ll^n,0|l(n) 



6''Cn ds 



(n-l+o) 



<(n-l + a)||0'Cn||(„+„)< 11^11(1) ||Cn||(„_i+„) 



Similarly, 



(n-l+a) 



< {n-l + a)B{a,n) \\0',\\^^^ lkn,o||(„) 



r(«)r(n) II 

r(n-l + a) ll''Hl{a)ll^«'Oli(n) 



and 







< 



(n) 



1(1) lpn,0||(„) • 



With these estimates we obtain 



// / n' 



6 Cnds 



1q I 9'^Znfids 



e'! / e'znflds 



r(a)r(n) II II II ,11 



< 



r(n + l + a)ll^oll«ll'^^ll(")"^'^'""W 
< B{l + a,n)\\e[\\^^^ lkn,oll(„)- 

(n+l+o) 



Finally, 



-l+a r(n) ||^,|| 



< u 



T{n + l + a) II ll(i)ll^"llW 



lly^llfn) 



and 



r(n + l + a) roii(i)ii^"'Oiiw 

Now we collect all the estimates from (|Sn|) through and obtain 



+ r(n + i + «) ('"^^"^ + ^"^^"^ + "^^^ + "^'1 
+ r(n + i + a) (ll^'lla) "^""w + "^"'Ollwj • 



(1) 
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This shows Cn+2 £ Fn+i+a- The above calculations and the bounds from Theorem |21 now 
imply the existence of a locally bounded function (j) : — > M"*" with (j){x) = 0{x) as 
X — > 0, such that with M = max{||6''||^]^j , ||^rll(Q)} ™d Q = (t){M) we have 

||Cn+2||(„+i+,) < (l + ^;^) IICn||(„_l+„) + Y^^^^YTa)^- 



Moreover, since 



lim n^_iM^ = 1 (58) 
n^oo r(n + /3) ^ ^ 

for each /3 > and ^'^r^J^ < 1 for each n e N and /3 > 1 (cf. [XKHtj . 6.1.46 and 6.1.47), 
the sequence (||Cn||(„+o_i))neN is bounded by the sequence (a„)„gN defined through 

2Q\ Q 



/ 2Q\ Q 



with a2 = ||C2||(i+a) < ll^rll(a)- '^n IS increasing, and so either (Q + l)a„ < Q for all n 
(then a„ < Q), or eventually ^^^^ < ""^^^"t'*'^ ; and then 

3(Q + i: 



. Q + l\ , 3(Q + 



This shows 

n/2 



where the infinite sum is bounded uniformly in a > inf > 0. The last inequality above 
follows by taking the logarithm of the product above and using ln(l + |3;|) < |x| in the 
resulting sum. Thus we obtain 

||Cn+2||(„+i+„) < max |</.(M), exp |^3(</.(M) + 1) j^^^ | , 

and the claim for ^„ is shown. 

Turning now to the ^n, (PH]) implies 

U = -KC'n-l - (^'oVn-l - e'rVn-l), (59) 

and gives 

Vn-i = - I {0[zn-i + e'^Cn-i) ds - y„_i(0) + y„_i,o(0). (60) 

JO 

The claim now follows in a very similar way as above from Propositions ^ and |21 □ 
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By the Lemma and Proposition^ a), 

11^ II ^11^ II r(n + l-a) ^ ( T{n + l-a) 

llUlloo ^ IIUII(„_l+a) fu+l-a - ^V^^l I tc-^-"r(n) 

Now ()58p implies 

(n - IV 

||Cn|L<c<^(M)^-^n-^+^ 

with 

V{n+l-a) i_„ 
c= sup — -TTr^TTT^'^ <oo, 
neN,aeJ<, r(n) 

and (|44jl is shown. The same reasoning apphes to ^n, showing (|43|1 and finishing the proof. 

□ 



4 General Hamiltonians 

Our main result Theorem ^ is formulated for Hamiltonians ((21) with constant eigenvalues 
satisfying Assumptions 1 or 2. In this section we show that these assumptions are satisfied 
for a large class of Hamiltonians after transformation to the natural time scale. 
Let us consider 

(leds - H{s))i^{s) = (61) 
for the traceless real-symmetric Hamiltonian 

~ ^ / Zis) Xis) \ ^ ( cose{s) sme{s) 
\X{s) -Z{s) J ^^ '\sme{s) -cos^(s) 

If X'^ + > 0, then for each Sf € M the transformation 



t{s) = 2 j ^/p^iu) dn (63) 

takes the equation ()61() with Hamiltonian (|62() into equation with Hamiltonian Q with 
9 = 9 o . Berry and Lim |BerLij found that under very general conditions on X and 
Z the singularities of 9' have the form of a first order pole plus lower order singularities. 
Then, as to be explained, by the Darboux' Principle Assumption 1 resp. Assumption 2 
are satisfied pointwise on the real line. More precisely, the n"^ derivative of ^' at t £ M 
behaves like T[n)r~^ as n ^ oo, where r is the distance from t to the nearest pole; the 
corrections have derivatives going like T{n — a)r~'^ as n ^ oo for some a > 0, and thus 
are in Fi_a,r- 

The task is now to make this discussion rigorous, and we start with giving a version 
of Darboux' Theorem. While this theorem in various forms is certainly well known jHel 
IBo| IDij , we were unable to find a statement in the precision and generality we need in the 
literature. The proof given in the Appendix uses Cauchy's formula and explicit integration 
near the singularities. This strategy was suggested to us by Vassili Gelfreich. 
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Theorem 5 (Darboux' Principle). Let f be analytic on Dn = {z ^ C : \z\ < R}, and 

assume that f is analytic also on dD^ except for finitely many points. Assume that there 
exists G N and (zj,aj,gj)j<:iy with the properties: 

(i) For each j, zj is one of the singularities of f on dDji; 

(ii) a,- EM\ {0,-1,-2,...}; 

(Hi) the function gj is analytic in a neighborhood Uj of z = 0; 

(iv) When zq is a singularity of f on dDn and Az^ := {j : Zj = zq}, then 

f{z)= ^{z-zor^'^gjiz-zo) on f] {Uj + zo) n Dr . (64) 

In words, Darboux' Theorem says that when / has finitely many algebraic convergence- 
limiting singularities, each of those contributes to the growth of the derivatives of / with 
a term of absolute value |(7(0)|r(n + Q)/i?"+" depending on the strength ^(0), order a 
and distance R of the singularity, and a phase depending on its strength and location. 
It is now clear that any function fulfilling the assumptions of Theorem [S] is in -Fa,_R({0}), 
where a = maxj aj. We have to prove this for the analytic continuation of 6' , and for this 
purpose put the following assumptions on the Hamiltonian H: 

Assumption XZ: X and Z are meromorphic in an open set U G C containing some 
point Sr on the real axis. X and Z fulfill 

p'^{s):=X^{s) + Z'^{s)>c>0 VsGMnC/. (66) 

By convention, we do not lift removable singularities of , so the critical points of 
consist of its zeros and the poles of X and Z . For s close to such a critical point sq of , 
we require 

X{s) = {s- sorfis - so) [l + {s- soTgxis - so)], 

(67) 

Z{s) = ±i(s - sorfis - so) [l + {s- so)"5z(s - so)-] 

where 0<nGN, mGZ and > —1; the functions f , gx o.n-d gz are analytic in a 

neighborhood of Sq, and K := \ f'^ifi){gx{^) ~ 9z{^))\ > 0. The set of critical points has no 
accumulation points in U . 

The class of X and Z fulfilling Assumption XZ is smaller than the universality class 
considered in |BerLij . It is not our ambition here to investigate just how large we can 
take our class while still giving a mathematically rigorous proof; but note that (|67() does 



Then 



/W(0) 



n! 
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Figure 1: The boundaries of the sets Cr for various situations; the black line is the boundary 
of Cr. In (a) and (b), = 1 + s^, and s^ = in (a) while s^ = 0.5 in (b). The branch 
cuts are on the imaginary axis. In (c) , = {1 + s^)^ and s^ = 0; examples (b) and (c) 
show why it is necessary to consider the connected components. 

contain the generic case of analytic X and Z and a simple zero of p^, i.e. m = and n = 1, 
along with many others. 

By the map s i-^ p{sY is analytic and invertible on f/riM, but may have critical 
points in the complex plane. For each such critical point sq for which a branch cut Bg^ of 
\/~f? is needed, we choose the branch cut such that it points away from Sr, and define Uq 
to be the connected component oi U \ \_]^^ Bgg containing s^. Then r as given in is 
well-defined and analytic in Uq. We define 

Cr{s,) := {seUo: |r(s) - T{sr)\ < r}, 

denote by Crfi the connected component of Cr containing s^ and put 

R = R{sr) := sup{r > : Cr,o(sr) C Uq}. 

Figure 1 shows some typical cases. 

We restrict the discussion to the case where Cr{sr) hits the boundary of Uq at critical 
points of as r grows to R. When U is large enough, this is the generic case, provided one 
has put the branch points Bgg in a sensible way. We distinguish two cases corresponding 
to our Assumptions 1 and 2 on 9': 

Assumption Rl: Let s^ £ ^ be a point where a Stokes line of t, i.e. a level line of 
-Re(r(s)), emanating from a critical point sq of p^ crosses the real axis. Assume that 
dCii{sy:) n dUo = {sqjJq}, which means that t(so) is closer to T(sr) than the r-image of 
any other critical point of p^ . 

Assumption R2: dCR^Sr) n OUq = {si, . . . s^}, where the Sj are critical points of p^. 

Theorem 6. Let H , p and 9 be defined as in \61\) with X and Z fulfilling Assumption XZ. 
Define 9 = 9 o r^^ with r given by il6'5|) . and to = t{sq). 

(i) If Assumption R2 holds, then 6' G -^i,_R({ir}), ^-e. Assumption 2 is fulfilled at t^. 
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(ii) Suppose that Assumption Rl holds, and that X and Y are given by ( |6'7| j at sq. Then 

±n 
2m+n+2 



with 7 = „ ^"10 there exists a closed interval I 3 s^ such that on I, 



where 9[ G Fa^R{I) for each a £ (0, 1) with a > 



2m+n 
2m+n+2 ' 



Proof. Let Assumption 2 hold; without loss we assume Sr = 0. Then r is analytic on Cr, 
while on dCu there are finitely many singularities. Let sq be a singularity, and let ^0 be 
the class of functions h which are analytic in a neighborhood of with /i(0) = 0. Then for 
s close to So, 

p\s) = 2K{s - so)''"+"(l + h^{s- so)) 
with hi £ Aq. Consequently 



T{s)-r{so) = 2/ (r - so)(2™+")/V2K(l + /ii(r - sq)) dr 



so 



4v2i^ , 2m+n + 2 

:{s - so)^^ (l + h^is - so)) (68) 



2m + n + 2 



with h2 £ Aq. 

Since by construction r has no critical points inside Cr, it is locally analytically in- 
vertible there. Since Dr := t{Cr) is the disc with radius R and center T(sr), global 
invertibility follows. Thus r is one-to-one from Cr, onto Dr with analytic inverse. 

We have 

e'(T(s)) = — V4 = — TT-^ arctan — s) = 5 (s), 

and taking s = r^^(t), t G D/j, shows that 6' is analytic on the circle Dr, fulfilling the 
first assumption of Darboux' theorem. Note in particular that X' Z — Z' X is non-singular 
on Dr by our convention of not lifting removable singularities of p^. For the behavior of 
9' near a singularity sq at the boundary of Dr, we revisit the calculation of Berry and 
Lim [BerLij . paying special attention to the error terms that arise. We have 

{X'Z - Z'X){s) = ±mK{s - so)""+"-'(l + h^{s - sq)), 

with /i3 € ^0 as above. Writing u = s — sq we now obtain 

±inKf72™+"-i(l + h^{a)) ±in(l + h^{a)) 



Hris)) 



2(2i^a2"^+")3/2(l + /ii(a))3/2 4y2Z^(2m+n+2)/2(l + hi{a))^/^ 

±m (l + /i3(a))(l + /i2(^)) 

(2m + n + 2)(t(s) - r(so)) (1 + /ii(a))3/2 
ibin 

T — r--. 7 — rrll + hA[a)), 

(2m + n + 2)(T(s) -r(so))^ ^ 
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with /14 E Aq, where we used ()68p m the second hne. It remains to find the form of 
near the singularity r(so). First note that since r is an integral, t{s) — r(so) = T{a). From 
(|(i8|) . t(c7) = Ka°'{l + h2{cr)) with the obvious K and a. The function 

a^cji^i/"(l + /i2((T))^/" 

is invertible in a neighborhood of ct = 0, and the inverse function h^{a) is an element of 
Aq. We then find 

2m+n+2 

/14(a) = /i4 o /i5 (^r(a)2/(2™+«+2)^ = ^ r(cT)ra+25.(r(a)) 

i=i 

with analytic functions gj. Putting things together and writing t = t(s), Iq = t{sq) we 
obtain 

, . / 2m+n+2 



^'">= (2^ + „ + 2)(t-t„) ('+ g (*-*„)-^s,«-«| (69) 



in a neighborhood of to- This has exactly the form (|64j) . and thus Darboux' Theorem 
shows (i). As for (ii), note that when Assumption Rl is fulfilled, by continuity of r still 
9Cr(s) n OUq = {so, So} for all s in a real neighborhood of Sr- All the above calculations 
only need information from the singularities, so they are valid without change, and the 
proof is finished. □ 

Example 1. (Landau Zener transitions): In the Landau-Zener model, X{s) = s, Z(s) = 
6 and consequently p{s) = VS'^ + s^. The critical points of p are at ±15, and 



R{0) = r(i(5) = 2 _ s2 ^ 

Jo 



s = . 



Moreover, at s = ±15 the functions X and Z have the form 1)6 7|1 with m = 0, n = 1, 
/ = ±16, gx = ±1/5 and gz = 0. Thus 

o'it) = ±^^^^ (1 + ± ^sf^^giit T 15) + itT ^S)^^^g2{t T i<5) 

where gi and 52 are analytic near ±15. Thus 



° ■ 3t-i6 3t + i6 

and O-c = 6' — O'q ^ Fx/3 5(1) for some / D {0}. In the simple situation at hand, it is easy 
to see that in fact 9[ G Fii^^^{I) for every finite interval /. When has more that one 
critical point on each side of the real axis, the situation is more involved and we refer to 
|BeTe2| . 
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Appendix 
Proof of Lemma [T] 

We start by converting (|36|) and (|37jl to a recursion for the Dn alone by plugging (|41jl into 
1)3 7|) . The result, in a somewhat expanded form, is 

|m/||2 71-1 

Dn = -^^^B{k,n-k)+ (70) 

k=l 

||/i/||2 n—1 / k—1 p. n— fc— 1 „ \ 

A;=2 \j=l i=l / 

n—1 fe— 1 „ n— fc— 1 „ 

+ ll^1la)E^(^'--^)Ef ^ T+ 

fc=2 j=l (=1 

n-2 

+ J^B(A:,n- A:)DfcZ)„_fc. (73) 

fc=2 

To show (|l2|) , we will of course proceed inductively. Direct calculation yields that (|l2|) is 
true up to n = 10 (even for M = 1). Let us now suppose that n S N is an even number 
and that (|42j) holds up to n — 2. We will show that (|42|) also holds for n and for this 
purpose treat each line of ((7n|l through ((75)) separately. We start with (f75|l. Using the 
induction hypothesis, we get 

fc=2 \ j=l / \ 1=1 

= (n-l)(n-2) Y.Bik-hn-k-l)Y^ ^ ^(1^ = (*0. 



fc=2 i=l Z=l 

|2 

1(1)' 



We sort this triple sum by powers p of ll^'ll?!-,, i.e. take p = j + I. The scheme is the 



following: 

p = 2: j = l ^ 1 = 1, A; = 2, 3,..., n-2. 

p = 3: j = l =^ 1 = 2, k = 2,3,...,n-A 

j = 2 1 = 1, fe = 4,5,...,n-2 



so for given p < n/2 (which is the highest power of that occurs) we have j running 

form 1 through p — 1, and for this j we have I = p — j and k = 2j,2j + 1, . . . ,n — 2[p — j). 
This gives 

n/2 n-2p+2j 

= (^^^m^T.Pti)^'{^ E Bik-l,n-k-l) 

p=2 \j=l k=2j 
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We will prove that ip{n,p) is bounded uniformly in n and p < n/2. We use the integral 
representation 

Bin,m)= [ n"-i(l-u)"-M'a (74) 







of the Beta function, and obtain 

,X p— 1 n—2p 



•^0 j=l fc=0 ^-^^ 



^i=i ^-^ ^ / \ fc=o 

The sum in the second bracket above is bounded by 2 uniformly on [0, 1], and thus ■ip{n,p) 
is shown to be bounded provided we are able to prove 

E (^) ^''"'(1 - ^)'^''"'^"'dn < 5. (75) 

To see (|75p . first note that the terms with j = 1 and j = p — I are equal to p/{2p — 3) 
and thus bounded by 2 since p > 2. For the remaining terms, we use u'^^~'^ < and 
(1 - ii)2(p-i)-2 < (1 - 'u)P-J.We then obtain 



Y,[ ] «''■"'(! - ^x)2(P-J)-2 < ^ r n^(l - n)P-^' = 1 



i=2 ^-^^ j=o ^-^ 

by the binomial theorem. This proves ((7SJ), and we obtain 

"/2|m/||2p^p ^||0'||2f,MP 

E3 ^ („-!)(„- 2) i:^^^2(;7rT)E^^ P<5) 

^ ' ^ ' p=2 ^ ' p=2 

since n > 12. We now turn to 1)72(1 . and start by proving 
Lemma 3. For m < n we have 

V ^ < 1 ^^'^ im-2l + l)\\e%^M^ 

^ j - m ^ (21- 1)11 ■ ^ ' 

j=i 1=1 ^ ' 
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Proof. We use the induction hypothesis to calculate 

j=l j=2 ^ Z=l /=1 j=2l •' ' 



The claim now follows from 



^ A;(A; - 1) j - 1 m m(j - 1) ' 

□ 



Using we get 



- ll^ll(i)A.(A:-l)(n-A;-l) (2j - 1),! 

(n-fc-2/)||g-||;|)M^ ^ 
^ ^ (2/-1)/! ■^*^^' 

Again we sort this by powers p of leaving the leading out. The scheme is 

p = 2: j = l ^ / = 1, A; = 3,4,...,n-3. 
p = 3 j = l ^ / = 2, A; = 3,4,...,n-5 
j = 2 ^ ^ = 1, A; = 5, 6, . . . , n — 3, 



and this time the general term is j = 1, . . . ,p — 1, I = p — j, k = 2j + 1, . . . ,n — 2{p — j) — 1. 

-iTT B(k,n-k) B(k-l,n-k-l) , , , - 

We use = ^„_i)(„_2) and obtam 



(*2) = , '\7^ E (79) 
(n-l)(n-2) ^ 



with 

p-l n-2p+2j-l 



*(".rt=E E (,J (2,-l)(2(p-j)-l) 



j=l k=2j+l 
p—l n— 2p— 1 



^0 \J / 
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The second sum above is obviously bounded by {n — 2p — 1) ^1^=1 ^^^^^ < 4(n — 2p — 1) 
on [0, 1/2], and since it is symmetric around u = 1/2, this bound is also valid on [0, 1]. On 
the other hand, the integral representation of the Beta function yields 

^^(P\ tx2i-i(l-n)2(p-j)-i J^(P\ B{2j,2{p-j) 



^ Vj7 (2J-1)(2(P-J)-1) ^ Viy (2j-l)(2(p-i)-l) 
^ (2.-l)W-2) g(')^^^^--^'^(^-^--^^)^ 



The last integral is bounded by 5 due to (f7S|) . and thus (j){n,p) < (2p^_?i)(^2p-2) — ^{p+i) ' 
Inserting in 1)79^ yields 

Turning to ()7ip. note first that for j < n — 1, we have 

n-l 

^ B{k,n-k) = B{j + l,n-j-l) + 
k=j+i 

+B{j + 2,n - j - 2) + ... + B{n - 1,1) = 
= B{k,n-k). 

k=l 

Moreover, 5(1, n- 1) = l/(n- 1) and B{k,n-k) < 2/((n - l)(n - 2)) for 2 < /c < n - 2, 
and thus 

n-l 

J2B{k,n-k) <4/(n-l). 

k=l 

By symmetry, 

n— 1 fe— 1 „ n— 2 „ n— 1 

jZB = 2||e'||^^)^i?(fe,n-fe)J^-^ = 2||0'||Jy j;-^ 5] i3(A;,n-A:)< 

fc=2 j=l j=l •' k=j+l 

I) D, ^ 8 ||g-||^,) (n - 2/) ll^^ll^l) M' 

n-l ^ i -(n-l)(n-2) ^ (2Z - 1)/! 

where the last inequality is (f77j) . Thus 

^ ' p=2 



< ■■ 1(1) 
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Finally 

1 n—l WO'W'^ 

® = 4 11^11 (1) E ^(^' ™ - ^) ^ ^rzf- 

fc=l 

Combining ((TBI), (IHOl), (IHU) and JHll) we arrive at 




Choosing M > 42, the proof of Lemma ^ is finished. 



Proof of Theorem [H 

Let zq be a singularity of / and write Uzq = Clj^AzQ^^i^' From (1^1]) it is clear that the 
function z i-^ X]jgyi^g(^ ~ -^o)""^ — ^o) is the unique analytic continuation of / from 
{Uzo + zo) n Dfj to {UzQ + Zq) \ Bzo, where B^q := {z e C : z = azo, a > 1} is the 
branch cut (if necessary). Moreover / is analytic on the closed set \ [Jj{Uj + zj), and 
therefore analytic in a neighborhood of that set. Putting this continuation together with 
the continuations near each singularity, we conclude that there exists 6 > such that the 
analytic continuation of / to D 11+25 \ [jj Bj exists and is bounded on (dDji+s) \ Uj ^j- 
Here we may choose 5 < 1 and sufficiently small to guarantee Ds G Uj and Dsr £ Uj for all 
j. Let r = dDji+s U U^i be the piecewise smooth path that encircles anticlockwise 
along the boundary of the disk with radius R + 5 and avoids the branch cuts B^^ by 
encircling clockwise the singularities at zj with a circle Cj of radius S. Then 



/W(0) 



nl 



27ri 7 27ri J 27ri J 



The first integral is easily estimated through 



2iri JaD- 



^ supi |^^_,,|/(Z)| ^j^^ , ^ 



For the contribution of the poles, first note that by (|64|) we may treat each j < N 
separately, even if they belong to the same pole. If aj S N, a straightforward computation 
shows that 

Res (M.) - (-l)".-iM2l!^ (1 + (1)) 

Noting that Cj is negatively oriented, we conclude that the contribution from poles is the 
one claimed in 

For the remaining terms let gj{z — Zj) = ^^Lo — zj)^. Then 

/ dz = > 6fc / ^-j-, dz = 

27ri J 27ri^ J 
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Cj I Zj 

where we substituted z = ZjC,. The remaining integral is shifted to the origin and can 
then be solved explicitly in terms of the hyper geometric function F, 



\k—ai 



J ^^1+1 dC = J (^^Ypl 

Cj/zj Cj/zj — 1 



jp, l + n,k-aj + l _^ 
' k- aj + 2 ' / r , 

5-Oi 



k — aj -\- 1 



^ ( 1 + n, A; - Q, ■ + 1 A 
k-a+2 



k — aj -\- 1 



In the second line above we used the power series expansion 



(1 + 0-+^ ,=0 

valid for |C| = 5 < 1. Note also that the branch cut was moved to the positive real axis 
through the two changes of variables. 

For k < max(a, 2 — a) we use the asymptotic expansion of the hypergeometric function 
(cf. HbStj, 15.7.2 ) 

P{^^k-~ + 2^ =(5nr (^--^+^)r(fc- a, + 2) + , 

which shows that from the first [max(a, 2 — a)\ + 1 terms in each sum only the A; = 
terms contribute to the leading order in ()65() with 
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Y{2 - aj) e-'™J sin((l - a^Ovr) _ n"^-^ e"^™ 



5,(0)- 



■3 



As to be shown, for k > max(a, 2 — a) we have 



Hence, these terms do no contribute to the leading order in ()65p . Note that the sum over 
k in (|55|) converges since -D^i^ i S C/j. 

To check (|84j) we use the integral representation of the hypergeometric function, 



- aj + 2 ' y r(l - aj + A;)r(l) Jo (1 + ^i) 



1 j-k—ai 



l + n,k- aj + 1 ^ \ _ r(2 - + /c) t 



\n+l 
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= ik-aj + l){ I dtj:^-^^^+ I dt- 



gk—aj J 
< (A;-a, + l)-— + 



1 1 

(1 + - Qfj + 1 



n(5 (1 + (5s)"+i 
[7?i (fe - ttj + 1) 1 _^ 1 _ / ^ ^ + 1 
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